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The recent Letter of Cowsik, Ratnam and Bhattacharjee claims that the best fit value 
of the velocity dispersion of the Galactic dark matter is ~ 600 km s -1 . If correct, this is a 
significant result that should lead researchers to design new experiments to detect very high 
energy weakly interacting massive particles. The interpretation of the Galactic microlensing 
data for massive compact halo objects would require substantial, far-reaching revision. 

The basis of the analysis of Cowsik et al. is an iterative scheme that solves the coupled 
Poisson and collisionless Boltzmann system of equations (Eqs. (3) and (4) of their Letter). 
There is no reason whatsoever to believe that the velocity distribution of the collisionless 
dark matter is Maxwellian. Rather, the distribution of velocities that supports the dark 
matter halo against gravitational collapse must be found by solving the self-gravitation 
equations. Almost all known dark matter halo models do not have simple Maxwellian 
velocity distributions ||. Dark matter halos with Maxwellians have a number of simple 
properties. The iso-density contours of the dark matter must coincide with the global 
equipotentials. For a flattened self-gravitating system, the contours of constant density are 
always flatter than the equipotentials ||. Only for spherical self-gravitating halos can the 
iso-density and equipotential surfaces exactly coincide. In the inner parts of the model 
described by Cowsik et al. , the halo feels the gravity of the disk and the spheroid and so is 
flattened. In the outer parts, the gravitational potential of the halo overwhelms that of the 
exponentially declining disk and the comparatively puny spheroid. As the halo feels only its 
own gravity field, the dark matter density distribution becomes round. The assumption of a 
Maxwellian velocity distribution now must inexorably drive the halo model to the spherical 
solution of the Poisson and collisionless Boltzmann equations with a Maxwellian. This 
always has the asymptotic behaviour of the isothermal sphere [[|. 

So, Cowsik et a/.'s numerical algorithm should yield a dark matter halo that deviates 
from the isothermal sphere only in the inner parts of the Galaxy. If the dark matter three- 
dimensional velocity dispersion is (v 2 )dm; t ne na l° must become indistinguishable at large 
radii from the isothermal sphere with rotation curve of amplitude \J\(v 2 )\^. Evidently, 
then, Fig. 1. of Cowsik et a/.'s original Letter Jl| cannot be the whole story. For example, 
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the model with (u 2 }dm ~ 600 km s™ 1 must match asymptotically onto an isothermal sphere 
with a rotation curve of amplitude ~ 490 km s -1 . But, Fig. 1 shows the rotation curve 
of this model is ~ 200kms _1 and apparently still slowly falling at a Galactocentric radius 
of ~ 20kpc. This can only be reconciled with our line of reasoning if the rotation curve 
rises very steeply at larger Galactocentric radii. That this indeed happens is confirmed 
by the astonishing rotation curve provided in Fig. 1 of Cowsik et a/.'s response || to this 
comment. This does have the right asymptotic behaviour - and betrays the physical origin of 
Cowsik et a/.'s high dark matter velocity dispersion! The model has an extensive envelope 
of dark matter at Galactocentric radii ~ IMpc. This is untenable, as it is in clear and 
unambiguous contradiction with the mass estimates of the Local Group provided by, for 
example, the timing argument J7|. To estimate the mass in Cowsik et a/.'s model, let us 
roughly approximate the rotation curve of Fig. 1 of as two piecewise continuous segments, 
one of amplitude ~ 250 km s -1 out to ~ 50kpc, the second of amplitude ~ 550 km s -1 from 
~ 50 kpc to 1 Mpc. This provides an estimate of the mass within 1 Mpc as 
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Peebles ]7J reckons the mass within the Local Group as ~ 4.3 x 1O 12 M , so Cowsik et al. 
violate this constraint by over an order of magnitude! It is the pressure required to balance 
the weight of the overlying layers of this phenomenal mass at large radii that causes the 
high dark matter velocity dispersion. Cowsik et al. argue that their result is robust for all 
pressure-supported haloes as the moment (w 2 )dm appears as the pressure term in the Jeans 
equations. This is untrue. The Jeans equations tell us that (v 2 )dm behaves roughly like the 
gravitational potential ifi. The deeper a potential well, the higher the velocity dispersion 
required to support the model. To achieve their anomalously high velocity dispersion, Cowsik 
et al. have simply taken a much deeper central potential than is reasonable or warranted! 
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